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We suggest a possible realization of a solid-state memory capacitive (memcapacitive) system. 
Our approach relies on the slow polarization rate of a medium between plates of a regular capacitor. 
To achieve this goal, we consider a multi-layer structure embedded in a capacitor. The multi- 
layer structure is formed by metallic layers separated by an insulator so that non-linear electronic 
transport (tunneling) between the layers can occur. The suggested memcapacitor shows hysteretic 
charge-voltage and capacitance- voltage curves, and both negative and diverging capacitance within 
certain ranges of the field. This proposal can be easily realized experimentally, and indicates the 
possibility of information storage in memcapacitive devices. 
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I. INTRODUCTION 

A recent experimental demonstration of a nanoscale 
memory-resistor (memristor for short) ^ has sparked nu- 
merous investigations in the area of materials and sys- 
tems that show history-dependence in their current- 
voltage characteristics^"—. This has also led to the re- 
cent proposal and theoretical investigation of memca- 
pacitors and meminductors, namely capacitors and in- 
ductors whose capacitance and inductance, respectively, 
depends on the past states through which the system has 
evolved*. Together with the memristor, the whole class 
of these memory-circuit elements promises new and unex- 
plored functionalities in electronics, and the combination 
of these memory devices with their "standard" counter- 
part may find application in disparate areas of science 
and technology, including the study of neuromorphic cir- 
cuits to simulate biological processes^. In this paper, we 
focus only on memcapacitors. 

Various systems are known to exhibit memcapacitive 
behavior including vanadium dioxide metamaterials*, 
nanoscale capacitors with interface traps or embedded 
nanocrystalsiSrJ^, and elastic capacitor o^^'^^ . As antici- 
pated in Ref. I3, memcapacitive effects may also accom- 
pany memristive effects in nanostructures, since in many 
of them the morphology of conducting regions changes 
in timei. Moreover, memcapacitors can be simulated by 
electronic circuits^. However, the number of experimen- 
tally identified memcapacitive systems is still very lim- 
ited. In addition, it would be of great importance to 
identify a mamcapacitive system that can be fabricated 
relatively easily, and is flexible enough to cover a wide 
range of capacitances. 

Here, we suggest a possible realization of such a sys- 
tem based on the slow polarization of a medium between 
a regular capacitor plates. There are several physical 
mechanisms that can potentially provide a slowly polar- 
izable medium. Examples include tunneling, activated 
drift of charged vacancies/impurities, slow ion penetra- 
tion through a membrane, etc. In this paper, we will 
consider the tunneling mechanism and discuss a solid- 



state memcapacitive system consisting of metal layers 
embedded into a parallel-plate capacitor as schemati- 
cally shown in Fig. [TJ In this realization, the inter- 
nal metal layers, together with the insulator material, 
form a "metamaterial" characterized by a long polariza- 
tion/depolarization time. The application of an external 
voltage to the capacitor leads to a charge redistribution 
between the embedded metal layers. The tunneling cur- 
rent between the layers depends almost exponentially on 
the applied voltage. This feature is important for the 
operation of our suggested device allowing for the "writ- 
ing" of information (in the form of medium polarization) 
with high-voltage pulses, and "holding" such information 
when low/zero voltages are applied. The resulting mem- 
capacitor exhibits not only hysteretic charge- volt age and 
capacitance-voltage curves but also both negative and 
diverging capacitance within certain ranges of the field. 
Due to its simplicity and unusual capacitance features 
we thus expect it to find use in both logic and memory 
applications. 

This paper is organized as follows. Sec. Ulldescribes the 
theoretical framework used to simulate solid-state mem- 
capacitors. Numerical simulations for the case of ac and 
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FIG. 1: General scheme of a solid-state memcapacitor. A 
metamaterial medium consisting of A'' metal layers embedded 
into an insulator is inserted between the plates of a "regular" 
capacitor. It is assumed that the electron transfer between 
external plates of the capacitor (with charge ±g) and internal 
metal layers (with charges Qk) is negligible. Therefore, the 
internal charges Qk can only be redistributed between the 
internal layers creating a medium polarization. 
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pulsed applied voltages are presented in Sec. IIIII and ITVI 
respectively. An equivalent circuit model is given in Sec 
IVl Finally, Sec. |VT] presents our conclusions. 



II. MODEL 

Quite generally, a circuit element with memory (mem- 
ristor, memcapacitor or meminductor) is defined by the 
relations^ 



y{t) = g{x,u,t)u{t) 
X = f{x,u,t) 



(1) 
(2) 



where u{t) and y{t) are any two input and output vari- 
ables (i.e., current, charge, voltage, or flux), g is a gener- 
alized response, a; is a set of n state variables describing 
the internal state of the system, and / is a continuous 
n-dimensional vector function. In this paper we will be 
concerned with only memcapacitors. Therefore, using 
Eqs. ([Tl[5|, we define a charge-controlled memcapacitive 
system by the equations 



Vcit) = C-'{x,q,t)qit) 
i = f{x,q,t) 



(3) 
(4) 



where q{t) is the charge on the capacitor at time t, Vc{t) 
is the applied voltage, and C is the memcapacitance, 
which depends on the state of the system and can vary 
in time. Below, we identify the internal state variables 
of the solid-state memcapacitor shown in Fig. [T] and 
demonstrate that such memcapacitor is indeed described 
by Eqs. (1511) . 

Let us then consider the system schematically shown 
in Fig. [T] consisting of N layers of metallic regions sepa- 
rated by an insulating material. This multi-layer system, 
of thickness (5, is embedded into a standard parallel-plate 
capacitor. For convenience, although this is not neces- 
sary, we assume that the insulating material in between 
the embedded metallic layers is the same as the one of 
the regular capacitor. We also assume that our memca- 
pacitor is designed in such a way that there is no electron 
exchange between the external plates and those embed- 
ded. This can be achieved by selecting an appropriate 
separation between the internal metal layers and the ca- 
pacitor plates and/or by using an insulating material be- 
tween the external capacitor plates and the embedded 
ones that generates a higher potential barrier. Therefore, 
in what follows, we assume that the tunneling only oc- 
curs between the internal metallic layers. Consequently, 
the total internal charge is always zero: X^feLi Qk{t) = 0, 
where Qk{t) is the charge on the internal metal layers at 
time t. 

The operation of the polarization-based memcapacitor 
relies on the redistribution of the internal charges Qk 
between the embedded layers caused by the electric field 
due to the charge q on the capacitor plates. Polarization 



of the metamaterial results in an internal electric field 
between the layers that is opposite to the electric field 
due to the plate charges. Therefore, for a given amount 
of plate charge, the internal charges tend to decrease the 
plate voltage or, equivalently, to increase the capacitance. 

Let us then calculate the internal charge distribution 
and consequent capacitance. A charged plane creates a 
uniform electric field in the direction perpendicular to the 
plane with a magnitude E ~ a/{2so£r), where a = q/S 
is the surface (of area S) charge density, Sq is the vacuum 
permittivity, and Sr is the relative dielectric constant of 
the insulating material. Using this expression, we find 
that, in the structure shown in Fig. ^ the external plate 
voltage is given by, 

Vc = 2dEg + dEi + [S-2Si]E2 + [S-2{5i+S2)]E3 + 



[S - 2 (5i 



5k-i)]Ek...-SE_ 



N, 



(5) 



where Eg = q/{2S£o£r) is the electric field due to the 
charge q at the external plate and Ek = Qk/{2SeoSr) is 
the electric field due to the charge Qk at the k-th embed- 
ded metal layer. Eq. ([5]) can also be written as follows 



Vc = 



Co 



•(A-2Afe_i) 



N 



2q 



Q 



(6) 



where A = S/d, A; = J2]=i ^j/d for i = 1, 2, iV - 1, 
Ao = 0, and Co = eoSrS/d is the capacitance of the 
system without internal metal layers. Therefore, the ca- 
pacitance of the whole structure is 



C 



2 Co 



2 -f Ef=i [A - 2A,_i] 



(7) 



From this equation it is already clear that, unlike a 
conventional capacitor, there may be instants of time in 
which the denominator of Eq. ([7]) is zero while the nu- 
merator is finite. This may happen when the internal 
metal layers screen completely the external field, despite 
the presence of a finite charge q on the external capac- 
itor plates. At these times one would then expect di- 
verging values of capacitance. In addition, Eq. ((7| does 
not enforce a positive capacitance. Indeed, as we will 
show in the examples below, at certain instants of time 
the internal metal layers may over-screen the external 
field, resulting in a negative capacitance. It is interest- 
ing to note that an hysteretic negative and diverging ca- 
pacitance has been found also in ionic memcapacitors, 
namely nanopore membranes in an ionic solution subject 
to external time-dependent perturbations^^. A negative 
capacitance has been experimentally observed in differ- 
ent solid-state systemsi^"— , but not accompanied by hys- 
teretic and diverging values of capacitance. 

In order to describe the dynamics of the internal 
charges Qk, let us define Vk — —Ek^k+i^k as the volt- 
age between k and k -\- 1 metal layers. The electric field 
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if eVk > U. In the above equations, h is the Planck con- 
stant. The tunnehng current density jk = Ik.k+i/S as 
a function of voltage drop Vk is depicted in Fig. [5] for 
two adjacent layers. A strong increase in current with 
increase of Vk is clearly observed. This is an important 
property for the operation of our suggested memcapaci- 
tor. 



III. AC VOLTAGE OPERATION 



FIG. 2: Tunneling current density as a function of voltage 
drop between two adjacent internal layers calculated using 
the following values of parameters: U = 0.5eV, 5k = 5nm, 
and m = me, with the electron mass. Inset: the energy 
level scheme. 



In this Section, we present results of numerical sim- 
ulations obtained for the case of a sinusoidal voltage 
V{t) = Vo sin(27r/<) of amplitude Vb and frequency / 
applied to a memcapacitor C connected in series with a 
resistor R. Such a circuit is described by the equation 



Ek^k+i between two neighboring layers is obtained by 
adding the electric fields due to charges at all layers and 
at the external metal plates: 



Ek,k+i — —2Eq — El 
-2q-{Qi + 



E, 



!■■• — Ek + Ek^ 

Qk) + {Qk+i - 



En - 
Qn) 



(8) 



The dynamics of the charge at a metal layer k is then 
determined by the currents flowing to and from that 
layer: 



dQk 
dt 



= 1, 



k-l,k 



Ik 



k,k+l, 



(9) 



where Ik,k+i is the tunneling electron current flowing 
from layer k to layer fc -I- 1 (note that for the top and 
bottom layers [k — \,N), there is only one term on the 
right hand side of Eq. (|9])). By considering the layer 
charges Qk as state variables, we immediately notice that 
Eq. © is similar to Eq. ^ . This demonstrates that the 
system under consideration is indeed a charge-controlled 
memcapacitive system. 

If U is the potential barrier height between two ad- 
jacent metal layers (see inset of Fig. [5]), the tunneling 
current induced by the voltage difference Vk between the 
two can be calculated using the following expressiona^i 



k,k+l 



Se 
2TrhSl 



U - 



eVk 



exp 



— —r^^ 



eVk 

[/ + ^ ) exp 



U 



eVk 
2 



(10) 



if eVk < U, and 



Ik,k + 1 



Se^V^ 
inhUSl 



exp 



eft. 14 



(12) 



where, for C, we use Eq. ([7]). Eq. (fT^ is solved numeri- 
cally together with Eq. © describing the internal charge 
dynamics. Below, we present results of numerical simu- 
lations for two different memcapacitor structures. In our 
simulations, a small value of i? = 10 is selected so that 
the voltage drop on the resistor is negligible. 



A. 2-layer memcapacitor 

The simplest system exhibiting polarization memory is 
a two-layer memcapacitor. In Fig. [3]we illustrate simula- 
tion results of such a device. Starting at the zero-charge 
state, the sinusoidal applied voltage induces electron tun- 
neling between the internal metal layers resulting in non- 
zero Qi and Q2- As it is shown in Fig. ISl^a), positive half- 
periods of V induce negative Qi and positive Q2 charges. 
In turn, these cause a screening electric field opposite to 
the electric field of plate charges. 

It is interesting that on the charge-voltage plot (Fig. 
|3l[c)), q{Vc) forms a non-pinched hysteresis loop. This is 
a quite unusual feature of a memory device, especially, in 
view of the fact that, at the present time, only memris- 
tive devices exhibiting pinched hysteresis loops have been 
observed experimentally. Physically, it is clear that when 
an internal polarization in the memcapacitor is present 
and the plate charge is zero {q = 0), the internal po- 
larization creates a non-zero voltage drop on the device 
Vc 7^ 0. Alternatively, this voltage drop can be compen- 
sated by plate charges, but then we get Vc — at q ^ 0. 
This explains why the curve does not pass through the 
(0,0) point. The corresponding capacitance hysteresis is 
shown in Fig. ^d). As expected, we find both nega- 
tive and diverging capacitance values. In the vicinity of 
Vc = 0, the capacitance changes from -l-oo to —00 at 
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FIG. 3: Simulation of two-layer memcapacitor with symmetrically positioned internal layers, (a) The charge on internal metallic 
layers and memcapacitor plates as a function of time t. (b) Voltage on memcapacitor, Vc, as a function of time t. Charge- 
voltage (c) and capacitance-voltage (d) plots, (e) Added/removed energy as a function of time t. These plots were obtained 
using the parameter values Vo = 7.5V, / = lOkHz, d = lOOnm, 5 — 66.6nm, S — 10~'*m^, £r — U — 0.33eV, R — IQ.. 




-7.5 -5.0 -2.5 0.0 2.5 5.0 7.5 

FIG. 4: Gharge-voltage plot at different applied voltage fre- 
quencies /. The decrease of the hysteresis at higher frequen- 
cies is a signature of memcapacitora^. The calculation and 
device parameters are as in Fig. |3l 

both sweep directions. In Fig[3]^d), these abrupt changes 
appear as two coinciding vertical hnes. 

Next, we consider the added/removed energy to/from 
the capacitor which is defined as 

i 

Uc= j Vc{T)I[T)dT. (13) 


From Fig. ISje) it is clear that the present solid-state 
memcapacitor operates as a dissipative device since the 
amount of added energy is on average larger than the 
amount of removed energy (resulting in positive values 
of Uc at all times). This occurs because the electron 
tunneling between internal layers is accompanied by en- 
ergy dissipated in thermalization processes due to the 



different electrochemical potential energies of metal lay- 
ers. In addition, in a real device we cannot exclude the 
phenomenon of local heating that would also contribute 
to dissipation of energji^^. 



In fact, the energy dissipation in electron transfer pro- 
cesses or the usual energy dissipation in dielectric ma- 
terials described by the imaginary part of the complex 
permittivity can be used as an alternative to our ap- 
proach. The heat released can drive a phase transition 
in a material between the plates, such as transition from 
crystalline to amorphous state and vice versa, accompa- 
nied by a change in material's permittivity. In this way, 
the memcapacitor acquires a long-term memory based di- 
rectly on the value of . The required material for this 
purpose can be similar to chalcogenide glass, which can 
be "switched" between two states, crystalline or amor- 
phous, with the application of heat, but possibly at a 
lower phase transition temperature. The device opera- 
tion protocol can then be based in applications of ac- 
modulated voltage pulses whose amplitude, duration or 
frequency can control the material's state. 



Finally, the frequency behavior of the charge-voltage 
hysteresis loop is shown in Fig. HI It is clearly seen that 
the hysteresis shrinks at higher frequencies. This is a 
typical behavior of memory devices^ related to the fact 
that at high frequencies the internal degrees of freedom 
of a memory device do not have enough time to respond 
to the external perturbation. Similarly, with increasing 
frequency, we have observed a decrease in capacitance 
hysteresis as well as in the rate of energy dissipation. 
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FIG. 5: Simulation of four-layer memcapacitor. (a) The charge on internal layers and memcapacitor plates as a function of 
time t. (b) Voltage on memcapacitor Vc as a function of time t. Charge-voltage (c) and capacitance-voltage (d) plots, (e) 
Added/removed energy as a function of time t. These plots were obtained using the same parameter values as in Fig. |31 The 
locations of internal layers are defined by parameters 5i = 0.0245, 62 = 33 — 0.488(5. 



B. Multi-layer memcapacitor 

The results obtained for multi-layer memcapacitors are 
similar to the two-layer memcapacitor case. The only 
interesting difference is that in the case of multi-layer 
memcapacitors, a richer internal charge dynamics can be 
observed. Charge dynamics is essentially determined by 
the distance between the layers and the potential energy 
barrier between adjacent layers. As a prototype of multi- 
layer memcapacitor, we consider a four-layer memcapac- 
itor. 

In our particular example two additional internal lay- 
ers (2-nd and 3-rd) were inserted between the two layers 
of the two-layer memcapacitor. The layers' positions are 
specified in the caption of Fig. [S] layer 2 is positioned 
close to layer 1, while layer 3 is placed in between layers 
2 and 4. As a result, the electron tunneling current can 
be high between 1-st and 2-nd layers which is reflected in 
high amplitudes of Qi and Q2- Concerning the charge on 
the third layer, it is essentially close to zero. We explain 
it by a tendency of maximum charge separation in these 
types of devices. Fig. [5] displays our simulation results 
for this case. One can also see that the capacitance hys- 
teresis and the added/removed energy (Fig. [S] (d) and 
(e), respectively) are similar to those of the two-layer 
memcapacitor as shown in Fig. [3l 

IV. RETRIEVING THE MEMCAPACITOR 
STATE 

At this point, an important question is how one can 
read the memcapacitor state. Here, we demonstrate that 
the most simple way to do it is by using a single voltage 
pulse. Small amplitude or short voltage pulses are not 



suitable for this purpose as they do not change the charge 
distribution on the internal layers considerably. There- 
fore, the measured value of capacitance, using small am- 
plitude or short pulses, would always be equal to Co, the 
capacitance in the absence of the internal metal layers. 
It is thus important to ensure that the applied pulse has 
a large enough amplitude and is sufficiently long. Con- 
sequently, the device state will be altered by the pulse. 
However, since the initial state is read by this measure- 
ment, this state can be, in principle, restored. 

Fig. [5] illustrates the reading of the internal layer po- 
larization in a two-layer memcapacitor. Quite generally, 
we start with a memcapacitor having a certain internal 
polarization and assume that at the initial moment of 
time the voltage drop on the memcapacitor is zero. Al- 
though this is not a steady state of the device (at i ^ oo 
and Vc = we expect q, Qk ^0), the memcapacitor can 
stay in such a state sufficiently long since the tunneling 
current is exponentially suppressed at low internal fields. 
Our idea is to apply a single voltage pulse to the device 
and monitor the amount of charge needed to make Vc 
close to the applied voltage V . This amount of charge 
Ag, as we show, depends on the device state. 

In our particular calculations, we consider two cases 
characterized by a different initial polarization. When 
Qi is negative at i = (Fig. [IK a)), the appfied volt- 
age pulse does not significantly change the device state 
(Qi, Q2) and the amount of charge Ag added to the 
plates is small (hence the capacitance is low). In the op- 
posite case, when Qi is positive aX t — {) (Fig. El^b)), 
much larger charge should be transferred to the capac- 
itor plates since "re-programming" of the internal state 
is required. Therefore, this measurement would monitor 
a higher value of capacitance. We emphasize that the 
measured capacitance value is determined by both the 
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FIG. 7: (a) Equivalent circuit model of two-layer memca- 
pacitor. Here, Ci and C2 are usual capacitors and i?" is 
a non-linear resistor. It is assumed that the internal layers 
are symmetrically embedded into the device. Otherwise, the 
capacitance values of all capacitors would be different, (b) 
Equivalent circuit model of A'^-layer memcapacitor. 



FIG. 6: The amount of charge Aq needed for the same change 
in Vc (by 7.5V in the present calculation) depends on the 
initial memcapacitor polarization, (a) and (b) represent the 
device dynamics at different initial conditions, while (c) is the 
voltage pulse profile. The memcapacitor parameters are as in 
Fig. El 



system state and by the voltage probe. For the same 
system state but different probes (for example, pulses of 
different polarity but of same magnitude) the measured 
value of capacitance may be different. 



V. EQUIVALENT CIRCUIT MODEL 

In this Section, we suggest an equivalent circuit model 
of the solid-state memcapacitor discussed above. This 
correspondence is very useful in circuit simulators, or in 
the realization of electronic circuits to reproduce memca- 
pacitive features. The main idea behind the circuit model 
is to represent each insulator spacing between adjacent 
metal layers (including capacitor plates) by a separate ca- 
pacitor and to use non-linear resistors in order to mimic 
the electron tunneling between the layers and energy loss 
in the thermalization processes accompanying the tunnel- 
ing. The top and bottom surfaces of each internal layer 
are considered as plates of adjacent circuit model's capac- 
itors while the metallic material of the layer itself plays 
the role of the conductor connecting these capacitors. 

Fig. [TJa) shows the equivalent circuit of a two-layer 
memcapacitor. The circuit is composed by three capac- 
itors connected in series. The capacitor in the middle 
is connected in parallel to a non-linear resistor (see Fig. 
[71[a)). The parameters of capacitors are determined by 
the usual parallel plate capacitor expression and related 



to the parameters of memcapacitor in the following way: 



Ci 

C2 



2 Co 



{d-S)/2 1 

EoEr S _ Cq 

6 " 'A' 



Next, we demonstrate the equivalence of the equations 
describing the circuit model and those of the memcapac- 
itor. The circuit shown in Fig. [Tl^a) is described by the 
following two equations 



Vc = 

dq 
dt 



2q & 
Ci ^ C2 
dQ^ 
dt 



(14) 
(15) 



where the first equation is for the total voltage drop and 
the second equation is the Kirchhoff's law for the cur- 
rents. Here, In stands for the current through the non- 
linear resistor which depends on the voltage drop on C2 
(equal to Q' /C2 with Q' the charge on capacitor C2). 

Solving Eqs. (|14I15|) for the total capacitance we ob- 
tain 



C 



Vc 



Ci C2 



Co 



2C2 



d{Q' 



1 + A'- 



dt 



-In 



(16) 
(17) 



In order to relate Eqs. (jl6ll7p to the memcapacitor 
equations from Sec. |lll we note that the sum of charges 
at the bottom plate of Ci (equal to —q) and at the top 
plate of C2 (equal to Q' , see Fig. [7[a)) is equal to Qi, 
that is the charge at the first layer of memcapacitor, i.e., 
Qi = Q' — q. One then immediately notices that Eq. 
(fT5)) and Eq. ([7]) (for the case of two layers) are exactly 
the same; Eq. p7)) and Eq. (jH]) are also the same. This 
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proofs the complete equivalence between memcapacitor's 
and circuit model's equations. 

The extension of the circuit model to an A^-layer mem- 
capacitor is straightforward: N — 1 different capacitors 
and non-linear resistors must be introduced between the 
external capacitors Ci. Fig. [Tljb) shows the equivalent 
circuit model of an A^-layer memcapacitor. Another in- 
teresting model extension can be achieved if a memristive 
dielectric material (e.g., VO2 films^) is used instead of the 
usual dielectric in the device. Then, in the equivalent cir- 
cuit model, the role of non-linear resistors will be played 
by memristors. 

VI. CONCLUSION 

In conclusion, we have suggested and analyzed a solid- 
state memcapacitor made of a multi-layer structure em- 
bedded in a capacitor. The resulting system has note- 
worthy features, such a frequency-dependent hysteresis 
under an ac-voltage, diverging and negative capacitance. 
These feature emerge due to the slow polarizability of the 
internal multi-layer structure as a consequence of electron 
tunneling between the internal metal layers. This gives 
rise to both complete screening of the field due to the ex- 
ternal metal plates - giving rise to diverging capacitances 
- and over-screening of the same field leading to negative 
capacitances. Clearly, when combined to an external cir- 
cuit with finite capacitance, the capacitance divergencies 
of our proposed solid-state capacitor would be cut off by 
the external circuit. Nonetheless, both the information 



stored in this memcapacitor in a continuous fashion (ana- 
log operation) and its negative capacitance in a certain 
range of the external field, may find useful applications 
in electronics. In particular, we have suggested a simple 
way to read the information stored in the capacitor by 
using appropriate single pulses. 

We have also shown the equivalent circuit model for 
this memcapacitor based on a combination of regular ca- 
pacitors and non-linear resistors. This analogy could be 
useful both in actual calculations with circuit simulators, 
or in experiments with electronic circuits to reproduce 
memcapacitive features. Finally, it is worth stressing 
that the memcapacitor presented in this work can be eas- 
ily fabricated by growing a metamaterial of several metal- 
lic layers separated by an insulator (e.g., an oxide) in be- 
tween the metallic plates of a regular capacitor. Since 
only "standard" materials (as opposed to ferroelectrics) 
are used, the suggested memcapacitor should also show 
high reliability. Several modifications of the suggested 
structure are possible, such as suggested in the text use 
of phase change or memristive materials as well as e.g. 
use of nano-size metal beads instead of metal layers. We 
thus hope this work will motive experimental research in 
this direction. 
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